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Abstract 

The Atlantic meridional overturning circulation (AMOC) is the northward flow of surface water to subpolar latitudes where 
deepwater is formed, balanced by southward abyssal flow and upwelling in the vicinity of the Southern Ocean. It is 
generally accepted that AMOC flow oscillates with a period of 60-80 years, creating a regular variation in North Atlantic sea 
surface temperature known as the Atlantic multidecadal oscillation (AMO). This article attempts to answer two questions: 
how is the AMOC driven and why does it oscillate? Using methods commonly employed by chemical engineers for 
analyzing processes involving flowing liquids, apparently not previously applied to trying to understand the AMOC, an 
equation is developed for AMOC flow as a function of the meridional density gradient or the corresponding temperature 
gradient. The equation is based on the similarity between the AMOC and an industrial thermosyphon loop cooler, which 
circulates a heat transfer liquid without using a mechanical pump. Extending this equation with an analogy between the 
flow of heat and electricity explains why the AMOC flow oscillates and what determines its period. Calculated values for 
AMOC flow and AMO oscillation period are in good agreement with measured values. 
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Introduction 

The Atlantic multidecadal oscillation (AMO) is the regular 
variation in North Adantic sea surface temperature (SST) with a 
period of 60-80 years first formally recognized (but not named) in 
1994 by Schlesinger and Ramankuttj [1] using spectrum analysis of 
temperature records from 1858-1992. Kerr [2] first used the term 
AMO in a 2000 review article, finking it to warming beginning 
around 1910 and peaking around 1940 followed by cooling ending 
in the mid-' 70s. Delworth and Mann [3] reported a temperature 
oscillation with a period of approximately 70 years using proxy- 
based reconstructions of surface temperatures during the past 330 
years and suggested fluctuations in the Atlantic meridional 
overturning circulation (AMOC) as the cause. Enfield et al. [4] in 
2001 reported a 65-80 year cycle in North Atlantic SST data for 
the years 1856-1999. Gray et al. [5] in 2004 created a tree-ring 
based reconstruction of the AMO going back to 1567 with a 60- 
100 year cycle and tied this to SST anomalies from Kaplan et al. 
[6] . In 2005 Sutton and Hodson [7] showed evidence for the AMO 
using an index based on annual mean SST observations and 
Polyakov et al. [8-9] reported a 50-80 year ocean surface 
temperature cycle that is exceptionally strong in the Arctic. 
Knudsen et al. [10] in 201 1 used climate proxy data to suggest that a 
55-70 year AMO has persisted over at least the last 8000 years. 
Knight et al. [11-12] and Wei et al. [13] claimed that the AMO is 
internally driven rather than externally forced. Hu et al. [14], 
Wang et al. [15], Kumar et al. [16] and Edwards et al. [17] noted 



the important influence of the AMO on regional and global 
climate. It has been observed that the increase in Northern 
Hemisphere surface temperature over the past century has not 
been mono tonic [18] and the AMO is believed to contribute to 
this variability [19-21]. 

It is generally accepted that the AMO is a result of regular 
variation in the AMOC [22-23]. The AMOC is sometimes 
referred to as the thermohaline circulation (THC), a term that wiU 
not be used here because of its imprecise meaning [24—25] . This 
circulation was first described in 1959-61 by Stommel and Arons 
[26-27], who proposed a global circulation comprising northward 
flow in the upper layer of the Atlantic Ocean ending in deepwater 
formation in the subpolar region, balanced by southward abyssal 
return flow to upwelling in the vicinity of the Southern Ocean. 
The AMOC is the primary means for transporting heat from the 
Atlantic equatorial region to northern latitudes [28]. 

Deepwater formation, the sinking of surface water due to 
cooling and ice formation increasing seawater salinity and density, 
was initially emphasized as the primary driver of the AMOC [29], 
but Kuhlbrodt et al. [30] in 2007 stated, "It is evident that surface 
buoyancy fluxes cannot provide the energy that is necessary to 
drive the AMOC." The problems in understanding how the 
AMOC is driven have been addressed by incorporating other 
processes in addition to deepwater formation, including ocean 
eddy fields and the influence of wind stress on overturning [31]. 
This has resulted in proposed AMOC circulation mechanisms of 
great complexity: driven by the surface buoyancy flux in the North 
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Figure 1 . Thermosyphon loop cooler. This apparatus is used in industrial applications for advective transport of heat from source to sink with a 
circulating liquid without employing a mechanical pump. Heat is added to the left leg and removed from the right leg. The density in the hot leg {pn) 
is lower than in the cold leg (pc), so the liquid level in the hot leg (Zh) must be higher than in the cold leg (Zc) to maintain the force balance between 
the two legs. The resulting difference in level (Zj,/) causes the flow to circulate. The apparatus is self-regulating because an increase in heat input 
increases the temperature of the hot liquid, which decreases its density, increases the difference in level between the two legs, increases the flow 
rate, and thus increases the rate of heat removal. The thermosyphon loop cooler behaves in basically the same way as the AMOC. 
doi:l 0.1 371/journal.pone.0100306.g001 



Atlantic and by energy from internal dilTusion in the upwelling 
branch and upper-ocean Southern Ocean to South Atlantic 
transport, all influenced by winds, basin geometry, bottom 
topography, advection, and small-scale processes [32]. 

Contrasting with this complexity is the idea that the AMOC is 
driven primarily by the meridional density gradient. This seems to 
have begun with Marotzke et al. in 1988 [33] who stated that 
circulation is driven by the density difference between equatorial 
and polar seawater. Rahmstorf in 1996 [34] claimed that the 
overturning circulation behaves as if driven by the Atlantic Ocean 
density gradient with flow having a linear dependency on density 
difference. Using an ocean general circulation model, Marotzke 
[35] indicated that the strength of the meridional overturning 
circulation is not much affected by wind forcing, except very near 
the surface. Straneo [36] mentioned that the loss of heat in 
northern latitudes creates meridional density gradients linked to 
the AMOC. Schewe and Levermann [37] indicated that 
meridional density differences ha\(; b("en assumed to be key to 
determining the overturning rate of the AMOC and listed several 
supporting references. Wang et al. [38] suggested that the AMOC 
flow depends on the density gradient between the subpolar and 
subtropical North Atlantic. 

An equation relating AMOC volumetric flow to density 
gradient is derived here that gives results in good agreement with 
measured AMOC flows. This equation is extended to relate 
AMOC heat flow to temperature gradient that is combined 
through an electrical analogy with AMOC thermal capacitance to 
give an equation for AJvIOC oscillation period that gives results in 
good agreement with observed AMO periods. 



Methods 

The methods used to understand how the meridional density 
gradient drives the AMOC wiU be familiar to chemical engineers 
because, in the chemical process industry, change in liquid density 
with temperature is often used as a means for inducing hquid 
circulation. These methods seem not to have been previously 
applied to analyzing how the AMOC is driven. It is best to begin 
with a description of an industrial apparatus that resembles the 
AMOC: the thermosyphon loop cooler. 

1. The Thermosyphon Loop Cooler and the AMOC 

The thermosyphon loop cooler is used in industrial applications 
where heat is to be transferred with a circulating liquid in 
situations where a mechanical pump should not be used because of 
the unacceptable consequences of pump failure or because of 
corrosive/hazardous properties of the circulating liquid. As shown 
by Fig. 1 , the thermosyphon loop cooler in its simplest form is a U- 
tube with a partially filled pipe connecting the tops of the two legs. 
Heat is added to the left leg and removed from the right leg, 
usually by way of integral heat exchangers (not shown). In some 
applications the heat source is an in situ exothermic chemical 
reaction and the circulating liquid is the reaction medium. Since 
the density in the hot leg (p//) is lower than in the cold leg {pc), the 
Uquid level in the hot leg {^) must be higher than in the cold leg 
i'Zc) to maintain the force balance between the two legs. The 
resulting difierence in level creates the driving force that 
causes the flow to circulate. The balance between the driving force 
and the opposing force created by friction determines the steady- 
state circulating flow rate, so the point of heat input "A" is located 
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as close to the bottom of the U-tube as practical and the point of 
heat output "B" as close to the top of the U-tube as practical to 
maximize the difference in level, and the cross-sectional area 
normal to the flow is made as large as practical to minimize the 
frictional force. The thermosyphon loop cooler is self-regulating 
because increasing heat addition increases hot leg temperature, 
decreases hquid density, increases difference in head {Zii^, and 
increases the flow of liquid (and the rate of heat transport) to the 
cold leg. 

The thermosyjohon loop cooler is like the AMOC in several 
respects: the hot leg corresponds to the tropical Atlantic with "A" 
being the bottom of the thermocline, the cold leg corresponds to 
the North Adantic deepwater formation region with "B" being the 
ocean surface, and the connecting partially-fiUed pipe corresponds 
to the North Adantic AMOC surface flow. But, in contrast to the 
thermo.syphon loop cooler, with the AMOC the bottom of the 
thermocline (a few hundred meters deep) is far above the bottom 
of the recirculation path (a few thousand meters deep). This is 
contrary to Sandstrom's 1908 argument [39] that thermally- 
induced circulation requires the heat source to be located below 
the cooling source at a depth close to the bottom of the 
recirculation path, but Jeffreys [40] in 192.^ rebutted that 
argument by showing how circulation can be maintained through 
temperature difiFerence irrespective of the relative heights of the 
heat source and sink, using the action of a thermosyphon loop as 
an illustration. 

Another difference between the thermosyphon loop and the 
AMOC is that the upwelling connecting the return flow from the 
cold leg to the hot leg is not directiy beneath the hot leg but many 
thousands of kilometers away in the Southern Ocean. This can 
happen because frictional resistance opposing return-flow through 
the deep ocean is relatively much smaller for the AMOC than for 
an industrial thermosyphon loop cooler since the ocean cross- 
sectional area normal to the flow is very large. The deep ocean acts 
as a reservoir, not a defined flow path, which is consistent with the 
lack of observational evidence for a direct connection between 
deepwater formation and abyssal return flow [41]. For the 
AMOC, the down-flow created by deepwater formation in 
subpolar northern latitudes is balanced by a corresponding up- 
flow at a location where the vertical temperature gradient is least 
supportive of density stratification: the Southern Ocean. The 
cross-sectional area normal to the up-flow is again very large, so 
the resistance to flow is very small. 

AH that is required to close the AMOC circulation loop is a 
means for conveying seawater from the Southern Ocean to the 
tropical Atlantic. This may be accomplished by the meridional 
component of the Adantic circumpolar current (ACC). The ACC 
is driven by wind shear from the strong westerly winds in the 
latitudes of the Southern Ocean, causing surface water to flow 
northward via Ekman transport resulting from the Coriolis force. 
The ACC flow rate is about 140 Sv [42-43] (one Sv= lO'' m^'s"'), 
much larger than the AMOC flow rate. Marshall and Radko [44] 
use residual-mean theory to describe how the ACC establishes 
meridional circulation and Rintoul et al. [45] discuss the 
relationship between the ACC and meridional circulation. 

Deepwater formation, although unable to drive the AMOC by 
itself, is a necessary condition for the existence of the AMOC since 
it is essential to the existence of abyssal return flow. The absence of 
significant overturning in the Pacific Ocean is consistent with the 
absence of extensive deepwater formation in that ocean. 

In contrast to the thermosyphon cooler, excess flow from any 
AMOC step can be diverted away from the succeeding step into 
the ocean proper, so the capacity for deepwater formation only 
needs to be equal to or greater than the rate at which the AMOC 



North Atlantic surface flow delivers seawater to the region of 
deepwater formation. This loose coupling among the three means 
that the slowest of the steps controls the rate of the overall process; 

so it is only necessary to understand the rate-controlMng step in 
order to understand the behavior of the AMOC. 

2. The Rate-controlling Step for the AMOC 

Figure 2 (a) shows how the density varies along the length of the 
AMOC flow path {Lpj, which begins in the tropical Adantic at the 
latitude corresponding to the beginning of the northward SST 
gradient slightly north of the equator. As shown by the 
temperature-depth curve of Fig. 2 (b), the upper temperature 
inflection defines the depth of the mixed layer {^j^a) and the mixed 
layer temperature {Tj^, and the lower temperature inflection 
defines the depth of the thermocline (Ztc) and the intermediate 
water temperature [Tm). The term "intermediate water" is used 
to identify the bottom of the thermocline and should not be 
confused with the same term when used to describe ocean water 
masses in other contexts. The mixed layer temperature is nearly 
constant with depth (within 0.5°C of the SST) while the 
thermocline temperature changes from Tml to Tjw in a nearly 
straight line with increasing depth. Below the bottom of the 
thermocline the slope of the temperature-depth curve is much 
smaller and gradually approaches deepwater temperature (2-3°C). 
The AMOC flow is driven by the difference in sea surface heights 
between the two ends of the flow path [Zdy)- 

The horizontal temperature and density gradients are at a 
maximum in the ocean mixed layer and decrease with increasing 
depth through the thermocline, becoming zero at the bottom of 
the thermochne. Moving north along the flow path, 7^/^ decreases 
while Thy remains constant. The point at which Tml equals Tjyy 
defines the end of the flow path for the present purpose. The 
extension of the horizontal density gradient beyond that point as 
the SST continues to fall below Tj^,-, eventually reaching the 
temp(;ratur(" of deepwater formation, is not relevant to the 
development of the equation for AMOC flow near the south 
end of the flow path. 

3. Equations for AMOC Flow and Period: Notation 

Named variables and constants that are used in multiple 
instances in the derivation of the equations for AMOC flow and 
oscillation period are listed here for easy reference. 

Linear dimensions (See Fig. 2). 

Zh Depth to the bottom of the thermocline at the flow path 

hot (south) end. 
Zc Depth to the reference level at the flow path cold 

(north) end. 
Zdif Difference between Zh and Zc- 
Zml Depth of the mixed layer at the ffow path hot end. 
Zjrc Depth of the thermocline at the flow path hot end. 
Zf Depth of flow path with vertically uniform density at the 

flow path hot end (calculated). 
Lp Length of the flow path. 
Wp Width of the flow path. 
Temperatures (See Fig. 2). 

Th Temperature (average) at the flow path hot end corre- 
sponding to 2^. 

Tiw Temperature at the flow path cold end and the 
thermocline bottom (intermediate water). 

Tjif Difference between Tu and Tiw 

Tml Temperature of the mixed layer corresponding to Zml- 

Physical properties. 

Ph Density (average) of seawater at the hot end of the flow 
path. 
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Figure 2. AMOC North Atlantic surface flow, (a) Darker color represents higher density. The hot (south) end of the flow path is defined 
(imprecisely) as the latitude where the SST northward gradient starts. The height of seawater above the bottom of the thermocline point (Zh) is the 
combined depth of the mixed layer (Zml) and the thermocline {Zjc)- The cold (north) end of the flow path is defined as the latitude where the mixed 
layer temperature becomes equal to the temperature at the bottom of the thermocline (7"/n/). At the cold end of the flow path the height of seawater 
(Zc) is defined as the vertical distance from the sea surface to the bottom of the thermocline at the hot end. The driving force for flow is the difference 
in sea surface heights between the two ends of the flow path {Z^jf). The vertical arrows represent the net heat flow out from the ocean surface due to 
radiative cooling, increasing with increasing latitude, (b) Typical tropical Atlantic temperature-depth profile representing the south end of the AlVIOC 
flow path. The temperature of the mixed layer at that point is nearly uniform down to the upper inflection point (Tml) and the temperature of the 
thermocline decreases linearly with depth down to the lower inflection point, marking the beginning of intermediate water (Tm). 
doi:1 0.1 371/journal.pone.01 00306.g002 



Pc Density of seawater at the cold end of the flow path. 
p Density of seawater averaged over the length of the flow 
path. 

Viscosity of seawater averaged over the length of the flow 
path. 

Cp Specific heat of seawater averaged over the length of the 
flow path. 

Pressure, force, velocity and flows. 

pjif Pressure difference equivalent to ^/f. 

F Force driving flow from d("nsity gradient. 

F' Force opposing flow from friction. 

V Average velocity across the flow path cross-section. 

V AMOC volumetric flow. 

qjw AMOC heat flow based on Tjiy reference temperature. 
Other terms. 
fjr Fanning friction factor. 

a Ratio of seawater meridional density gradient to merid- 
ional temperature gradient. 

JH Change in AMOC heat content per unit volume over the 
length of the flow path. 



G Thermal conductance (advective) of the AMOC volumet- 
ric flow. 

C Thermal capacitance of the AMOC heat storage layer. 

R Coefficient in the exponential excitation function. 

s Frequency in the exponential excitation function. 

t Time in the exponential excitation function. 

T Period of oscillation. 

4. Equation for AMOC Flow 

The steady-state AMOC flow must balance the driving force (a 
function of meridional density gradient) against the opposing 
frictional force (a function of velocity). The seawater densities at 

the hot and cold ends of the flow path (pjj and pc) ^"<i the length 
of the flow path [Ljj) determine the density gradient: 
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The vertical force balance requires that 

PcZc = PhZh or Zc = Zh{phIPc) so 
Zc = PrZhI {Ph + VpLj-) 

The difference between Zh and Zc is 

Zdif = Zh — Zc 

= Zh-PhZh/(Ph + ^pLf) (1) 
= Zh[1-Ph/{Ph + '^pLf)] 

Converting /^/j/^to pressure difference, using p for the average 
density over the length of the flow path, gives 

Pdif = PgZnll -PhUPh + VpZ-f )] 



The force from the meridional density gradient that drives 
AMOC flow is the product of the pressure difference and the 
cross-sectional area normal to the direction of flow. The width of 
the flow path (Wpj is bounded by the ocean width, but there is no 
physical boundary separating moving from stationary liquid that 
can be used to define the depth of the flow path. The development 
of the equation for the force opposing flow to follow assumes 
uniform flow over the cross-sectional area (vertically and 
horizontally), which is different from the actual situation for the 
AMOC where the horizontal density gradient driving flow is a 
maximum within the mixed layer [Zml) and declines to zero with 
increasing depth through the thermocline [Zt<:})- To correct for this 
it is necessary to calculate a depth [Zf) for a flow path ha\'ing 
vertically uniform density that exerts the same horizontal force- as 
the actual mixed layer plus thermocline. Rearranging Eq. 1 and 
substituting mixed layer (pmli and intermediate water densities 
ipiw) for the corresponding hot and cold densities gives: 

Zp = Zdif /(I - PmlIPiw) 

The value of Zf will be smaller than that of Zh (where 
Zh~ Zml+Ztc) depending on the relative depths of the mixed 
layer and thermocline {Zml and Ztc)- The force driving AMOC 
flow is then: 

F = ( WFZF)pdif = pgWpZp^ll - PhUph + V/9Lf )] (2) 

The frictional force opposing flow is proportional to the area of 
the wetted surface separating the moving and stationary liquid and 
the kinetic energy per unit volume of the liquid (Bird, Stewart and 
Lightfoot [46]). For a flow path of width Wp and length Lp the 
area of the wetted surface is WpLp. The kinetic energy term is 'A 
pv^. The proportionality factor relating these terms to the actual 
force is the Fanning friction factor, fp. 

F'=fpiWpLp){'/2pv^) or fp = 2F'/{WpLppv^) (3) 



This equation is not a law of fluid dynamics; its only purpose is 
to define the friction factor. An alternative to the use of Eq. 2 for 
defining the friction factor is the Darcy-Weisbach equation, a 
phenomenological equation having the same form. The result is 
the Darcy (or Moody) friction factor (used by mechanical and civil 
engineers), which is four times larger than the Fanning friction 
factor (used by chemical engineers), a frequent source of confusion. 

For laminar flow, the Hagen-Poiseuille equation, which is an 
exact solution of the Navier-Stokes equations, relates the friction 
factor to the Reynolds number. The AMOC is in laminar flow, 
even though the ocean is usually referred to as being turbulent. It 
would be better described as agitated. With turbulent flow the velocity 
fluctuates chaotically with time at every point throughout the 
volume of liquid and the energy causing this is derived directly 
from the flowing liquid. With agitation the energy creating chaotic 
velocity fluctuations is injected into the liquid from an external 
source (a turbine or propeller for industrial processes, wind shear 
for oceans) and the effects of this (-nergy injection diminish in 
proportion to the distance from the energy source. Most of the 
energy injected into the ocean by wind shear is dissipated within 
the Ekman layer, a depth of only 10—20 m. From the Hagen- 
PoiseuHle equation: 

fF = l6/Re=l6n/{Dhvp) 

In the above equation Z)/, is the hydraulic diameter, a parameter 
that allows the equation to be applied to flow c:ross-sections that are 
not circular. The definition of hydraulic diameter is D/,— 4A^^/P^^ 
where A^^ is the cross-sectional area normal to the direction of flow 
and -P„„ is the perimeter of the wetted surface. 

Since A,,^ Wp Zf and, for Wp»Zf, Pws^ Wp, dien D^^ 4Zp 
and the expression for the friction factor becomes 

fp = \6ti/{AZpvp) =Aiil{Zpvp) (4) 

Combining Eqs. 3 and 4 gives the equation for the opposing 
force: 

F' = 2h{WpLf/Zp)v (5) 

The equation for volumetric flow can now be obtained from the 
balance between driving and opposing forces, F'=F. Setting Eqs. 
2 and 5 equal to each other and solving for velocity gives: 

v= [pgZFy{2,x)][\-pH/{PH + '^pLF)]/Lp (6) 

Multiplying velocity by the flow cross-sectional area (Wp Zf) 
gives the equation for AMOC volumetric flow: 

^= 1/2 [pgWpZp'/fi] [I -pjj/{pH + YpLp]/Lp 

Since P}j»Vp Lp and Pn^p, the equation for volumetric flow 
simplifies to 

^={i/2gWFZpyn)Wp (7) 
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5. Equation for AMOC Oscillation Period 

The equation for AMOC period derived in this section depends 
on the analogy between heat and electrical flows. The first step in 
developing the equation is to com-ert the expression for AMOC 
volumetric flow as a function of density gradient to an expression 
for AMOC heat flow as a function of temperature gradient. 
Assuming an approximately linear inverse relationship between 
temperature and seawater density over the temperature range of 
interest gives Eq. 8, where a is the ratio of density gradient to 
temperature gradient, Vp/VT, 

S ( y^ga WfZp'^Ih) (Vr) (8) 

Equation 8 is multiplied by the change in heat content per unit 
volume, AH= p Cp [Tjj -Tm!), to obtain AMOC heat flow. 

Substituting the difierence in temperature between the two ends of 
the flow path divided by the length of the flow path {Tj,y-/Lp) for 
VT, gives Eq. 9, where G is the advective thermal conductance of 
the AMOC. Since qnv— 'P '^U, values for qjiy can be obtained 
directly from measured values for !f . 

qiw = ( '/2ga WfZ/AH/h) {TsflLp) = GT^if (9a) 

G = 'l2ga WpZ/AH/i^iU) = qiw/Tdif (9b) 

Notice that heat flow equals conductance times the temperature 
difference, which is Ohm's law. For .steady-state AMOC heat flow 
this is a trivial result, but it suggests the possibility of an electrical 
analogy for describing the dynamic behavior of the AMOC flow. 
This is influenced not just by the meridional temperature 
difference (voltage) and opposing frictional force (resistance) but 
also by time-dependent variation in the heat stored by the ocean 
upper layer (capacitance). The thermal capacitance of the AMOC 
is equal to the product of seawater density, specific heat and the 
volume defined by the width and length of the flow path and the 
combined depth of the mixed layer and thermoclrne: 

C = pCpWfLpZh (10) 

The electrical analog to the AMOC is a resistor and capacitor 
connected in series so the current (heat) flow through the resistor 
equals that through the capacitor. At time zero (south end of the 
flow path) the capacitor is fuUy charged. Increasing time 
(movement northward along the AMOC flow path) corresponds 
to decreasing charge in the capacitor (loss of heat by radiative 
cooling). When the capacitor is fully discharged (north end of the 
flow path) the total current (heat) flow through the capacitor, 
which is capacitance times the rate of change in voltage 
(temperature) with time, is equal to the total current (heat) flow 
through the resistor: 

C{dT,if/dt) = GT,if (11) 

Substituting the exponential excitation fimction for the 

relationship between voltage (temperature) and time, Tj,y= K e' 
where t is time and s is the complex-frequency variable (called this 
because it can assume real, imaginary or complex values and has 



dimensions of reciprocal time) [47], and solving the resulting 
differential equation gives 

Ke"{sc-G) = Q 

This equation is true for either K= 0 or sC= G. The first case 
corresponds to an initially uncharged state. The second case is 

s=G/C 

This is the natural frequency of the circuit (radians per second), the 
unforced frequency at which it will oscillate in response to a 
disturbance. It will also be the natural frequency of the AMOC. 
The corresponding natural period (seconds) is 

T=2tiC/G (12) 

Equations 7 and 12 together enabk- calculation of AMOC flow 
and oscillation period. What remains is to test the results obtained 
from these equations against observed values. 

Results and Discussion 

1. Independent Variables 

The independent variables that determine AMOC flow in Eq. 7 
are: the density gradient (Vp), flow path width (114) and depth to 
the bottom of the thermoclrne [Zh- Zml+Ztc)- The density 
gradient is a function of the difference in densities between the two 
ends of the flow path and the length of the flow path (Lpj. The 
density at the south end of the flow path {pfj} is a function of Tml, 
TjWi Zmu imd Zrc- This is calculated by assuming a linear 
temperature profile for the thermocline between Tml ^it depth Zml 
and TjiyaX depth ^and determining seawater densities (constant 
salinity) at five intermediate temperature points using a six-term 
polynomial approximation of the UNESCO equation of state for 
seawater. This is combined with the values of ZiiL and Zrc to 
calculate a depth-weighted hot-end average density (p/j). The 
density at the cold end of the flow path [pc) is determined from the 
temperature at the bottom of the thermocline {Tjiyj, which applies 
over the fuU depth of the flow path at that location. 

So the calculated values for AMOC flow CF) and oscillation 
period (T) are determined by six independent variables: 
temperatures of th(' mixc'd layer and intermediate water at the 
base of the thermoclrne [T^l and Tjj^, depth of the mixed layer 
{ZmlI, depth to the bottom of the thermocline (,^), and the length 
and width of the flow path [Lp and Wp). Values (or ranges) for the 
independent variables are: 

Tml=2& to 32°C 

r^= 10°C, from Fig. 2 (b) 

= 60 to 140 m, from Chu & Fan [48] 

Zh= 300 to 400 m, from Fig. 2 (b) 

The values for the preceding variables are based on measure- 
ments and are reasonably certain. But the values for flow path 
length and width are fuzzy since the actual flow path is very 
complex and prone to wandering. About all that can be demanded 
is that the values be consistent with measured distances between 
the pertinent SST isotherms and the physical width of the North 
Atlantic. With values set for the first four independent variables 
{Tml, Tjiy, Zml> and Zh) along with the desired values for AMOC 
flow and period (V and T), then finding a value for either flow 
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path length [Lpj or width (Wpj consistent with that result 
determines the required value for the other. 

A reasonable location for the hot (south) end of the AMOC path 
is just north of the equator where the northward temperature 
gradient begins, around 10 to 20°N. This is not far from the 
26.5°N location where AMOC flow is being monitored. The 
equator itself is not a suitable location for the hot end of the flow 
path because the meridional temperature gradient at that point is 
essentially zero. A first estimate for path length [Lj^ is the distance 
from the hot end of the flow path to the location of the SST 
isotherm corresponding to the temperature at the base of the hot- 
end thermocline (10°C), around 50 to 60°N, based on SST profiles 
(Space Science and Engineering Center, University of Wisconsin- 
Madison, http://www.ssec.wisc.edu/data/sst/). Considering sea- 
sonal variation, the latitudinal wandering of the isotherms, and the 
sinuous nature of the AMOC flow, only very rough estimates of 
path length can be made. With the four independent variables 
other than path length and width fixed and the desired values for 
AMOC flow and period chosen, selecting a value for path length 
also fixes the value for path width [Wfj. This is illustrated by Fig. 3, 
showing values for flow path length and width consistent with 
values for AMOC flow (tf') and period [T] of 17 Sv and 65 years 
as a function of mixed layer depth {Zml) for different values of 
depth to the bottom of the thermochne [ZhI- The path widths 
required for an AMOC flow of 1 7 SV and period of 65 years 
change in the opposite direction from the required path length 
with increasing mixed layer depth. 

2. Reality Checks 

As a first reality check, Eq. 1 is used to calculate the difiFerence in 
head [Zdijl between the two ends of the flow path. Variation in 
ocean surface height due to wind-driven (geostrophic) flows is less 
than two meters by satellite altimetry, so ^ifvcmst be smaller than 
this measured variation in ocean surface height but not so small as 
to be lost in the noise. The calculated values for range between 
0.4 and 0.7 m over the chosen range of values for the independent 
variables, a reasonable result. 



An independent reality check is to confirm that the power 
available from the AMOC as a heat engine is consistent with the 
pumping power required by the AMOC volumetric flow. All heat 
engines, natural or engineered, operate by adding heat from a 
source (solar heat gain) to a working substance (seawater) that is passed 
through a working body (the AMOC) to a heat sink (radiative heat 
loss). The working body converts part of this heat flow to power by 
exploiting certain properties of the working substance (change in 
seawater volume with temperature). The Malone engine [49-50], 
developed around 1925, created rotary motion by exploiting the 
change in water volume with temperature (it never attained 
commercial success). The pumping power required by the AMOC 
is the product of pressure difiFerence and volumetric flow. A head 
of one-meter equals a pressure of about 10 kPa, so for an AMOC 
flow of 18 Sv the power required is about 0.2 TW, a tiny fraction 
of the AMOC heat flow of about 1200 TW. Therefore, die 
pumping power required by the AMOC flow can easily be 
extracted from its heat flow, even if it is a very inefficient heat 
engine. 

There are two assumptions in the derivation of the AMOC flow 
equation that can be confirmed here as part of reality checking. 
The first is that the AMOC flow is laminar; the velocity («) 
calculated from Eq. 6 gives a Reynolds numbers between 10 and 
40, far down in the laminar range. The second is that the ratio of 
density gradient to temperature gradient (ot) is approximately 
constant; the calculated value of a over a temperature range of 26— 
32°C is between —0.23 and —0.26, validating this assumption. 

3. Calculated vs. Observed Values for AMOC Flow and 
Oscillation Period 

The most important test of the derived equations is the 
comparison between calculated and observed values for AMOC 
volumetric flow and oscillation period. Figure 4 (upper) shows 
calculated AMOC flow (If) as a function of the mixed layer 
temperature {TmlI for mixed layer depths [Zml) of 60, 80 and 
100 m. Flow path length {Lp) and width {Wpj are both set to 4 xio'' 
m and mixed layer temperature (Tml) is set to 10°C. Recent 
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Figure 3. Values for flow path length and width to give a calculated AMOC flow of 17 Sv and oscillation period of 65 years. Mixed 
layer and intermediate water temperatures (T^l and Tiw) are set at 28 and IO C. IVlixed layer depth (Z^/.) varies from 60 to 140 m and depth to the 
bottom of the thermocline (Zh) is set at 300, 350, and 400 m (blue, red and black lines). Dashed lines are flow path length (Lf); solid lines are flow path 
width (Wf). Increasing mixed layer depth corresponds to increasing flow path length combined with decreasing flow path width to give a calculated 
AMOC flow (V) and oscillation period (J) of 17 Sv and 65 years. There is no practical way of accurately determining actual values for either flow path 
length or width, but the values shown are reasonably consistent with distances between the North Atlantic SST isotherms corresponding to the 
mixed layer and intermediate water temperatures, and with ocean widths. 
doi:1 0.1 371/journal.pone.01 00306.g003 
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Mixed layer temperature {T„J - °C 

Figure 4. AMOC flow and period vs. mixed layer temperature. 

The curves show calculated AMOC volumetric flow CF) and oscillation 
period (7) for mixed layer temperatures at the south end of the flow 
path (Tml) ranging from 26 to 32'C with mixed layer depths (Z^t) of 60, 
80, and 100 m (blue, red and black lines). The curves relating flow to 
mixed layer temperature show the values for density gradients (Vp) 
calculated from the equations. The depth to the bottom of the 
thermocline (Z^) is set at 350 m. The length of the flow path (L^) is set at 
4000 km, roughly the annual average distance between the 28°C SST 
isotherm and the 10°C SST isotherm. The width of the flow path (Wp) is 
set to the same value, consistent with Fig. 3. The shaded areas on the 
figures indicate the approximate range of observed values for AMOC 
flow (1 2-22 Sv) and oscillation period (60-80 years). Values for flow and 
oscillation period cannot be calculated independently; fixing one 
necessarily fixes the other. 
doi:1 0.1 371 /journal.pone.01 00306.g004 

published values for flow from the RAPID AMOC observing 
system [51] located at 26.5°N are in the range of 12 to 22 Sv, 
consistent with calculated values. The flow rate is a function of the 
meridional density gradient (Vp), values for which are shown on 
the upper figure. The lower part of the figure shows the 
corresponding calculated values for the AMOC oscillation period 
(T). Observed values for AMO periods are in the range of 60-80 
years, also consistent with calculated values. 

Figure 5 shows calculated AMOC flow {1') and period (T) as a 
fimction of depth to the bottom of the thermocline {^„t) at the hot 
end of the flow path for mixed layer temperatures {T^l} of 28, 30 
and 32°C. Flow path length {Lj,j and width {Wp) are both set to 4 
xlO^ m and mixed layer depth is set to 80 m. 

4. Nature of the Equations 

The equations are derived entirely from first principles: there is 
no empirical or phenomenological content nor any ad hoc factor 
that can be adjusted to improve the agreement between calculated 
and measured results. Once values for the equation independent 
variables are set, so are both AMOC flow and oscillation period 
since the dimensions of the flow path (depth, width and length) are 
common to the calculations of both these quantities. Judicious 
selection of values for the equation independent variables that lie 
comfortably within their range of uncertainty can result in 




300 320 340 360 380 400 

Depth to the bottom of the thermocline 
(Z„) - m 

Figure 5. AMOC flow and period vs. depth to the bottom of the 
thermocline. The curves show calculated volumetric flow {f) and 
oscillation period (T) for depths to the bottom of the thermocline at the 
south end of the flow path (Zh = Zml+Z.tc) ranging from 300 to 400 m 
for mixed layer temperatures (J^d of 28, 30, and 32 C (blue, red and 
black lines). The length and width of the flow path (Lp and Wp) are the 
same as for Fig. 4. The sensitivity of the calculated values for AMOC flow 
and period to the depth to the bottom of the thermocline is high 
because the depth of the flow path Zp used in the calculation of AMOC 
flow is related to the depth to the bottom of the thermocline, ranging 
between 0.5 and 0.7 of Z^ depending on the relative depths of the 
mixed layer and thermocline. This term appears as Z/ in Eq. 7. The 
shaded areas on the figures indicate the approximate range of 
observed values for AMOC flow and oscillation period. 
doi:10.1 371/journal.pone.01 00306.g005 

remarkably good agreement between calculated and observed 
values for AMOC flow and oscillation period. Although supportive 
of the basic soundness of the equations, this should not be 
mistaken for an ability to accurately predict the behavior of the 
AMOC in the absence of observational data uncertainty that can 
be used for calibration of the independent variables. 

Although the AMOC is generally seen as a very complex 
phenomenon, the density-driven flow described by the derived 
equations is not complex. The main contribution to the apparent 
AMOC complexity must be the superimposed wind-driven 
geostrophic flows. These are comparable in magnitude to AMOC 
density-driven flow, but evidently their net contribution to 
meridional flow is small. There are also secondary deep flows 
created by viscous drag, causing the actual velocity profile to 
extend beyond the bottom of the flow path depth used in the flow 
calculation and beyond the thermocline as well. Because of total 
energy constraints these secondary flows are not likely to have 
much effect on total AMOC volumetric flow. 

The methods used here for deriving the equations are widely 
used by chemical engineers, but seem not to have been previously 
applied to the AMOC. Despite the simplifying assumptions 
implicit in the equations, especially the assumed confinement of 
the AMOC to the defined flow path and the neglect of horizontal 
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or vertical velocity gradients, they do a very good job of predicting 
the behavior of the AMOC/AMO. But they may not have much 
applicability to other ocean flows since the AMOC may be 
uniquely simple in the way it is driven. 

Conclusions 

1 . The AMOC is Driven Directly by the Meridional Density 
Gradient and Indirectly by the Corresponding 
Temperature Gradient 

The AMOC flow calculated from equations based on the 
meridional density gradient driving force agrees well with 
measured flow. 

2. AMOC Flow Oscillates at its Natural Frequency to 
Create the AMO, which is an Oscillation that is a 
Permanent Feature of the Earth's Climate System, not a 

Stochastic Variability 

The natural oscillation frequency calculatc'd from the advcctive 
thermal conductance of the AMOC and ocean upper layer 
thermal capacitance agrees well with observed AMO frequency. 
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